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I.  INTRODUCTION. 


The  purpose  of  this  report  is  to  validate  the  Sommerfeld  subroutine  in  our 
CDC-6500/6600  NEC  program,  and  to  present  the  results  for  the  standard  HF 
first-resonant  dipole  antenna  in  graphic  form  for  analysis,  reference,  and 
subsequent  comparisons.^. 

Our  original  AMP  (Antenna  Modeling  Program)  used  the  Fresnel  RCM  (Reflec¬ 
tion  Coefficient  Method)  to  obtain  antenna- to-ground  mutual  coupling  imped¬ 
ance,  and  it  has  been  pointed  out  that  a  Sommerfeld  method  must  be  used  when 
the  height,  H^,  of  any  part  of  the  dipole  over  earth  is  less  than1  '2 


0.70 


wavelength 


where. 


e  =  e  -  j 

r 


1.79751  x  104  T 
fMHz 


numeric 


*  earth's  relative  dielectric  permittivity;  numeric 
T  =  earth's  conductivity;  mhos/meter 

Using  the  HF  frequency  range  together  with  earth's  electrical  properties 
listed  in  Table  I,  equation  1  is  plotted  on  Figure  1. 

The  results  plotted  on  the  next  50  graphs  are  mutual  impedance,  R21  and 
X2 i ,  solutions.  These  solutions  were  obtained  by  using  NEC  to  compare 
the  antenna's  self  (free-space)  impedance,  Zj  1  ,  with  the  antenna's  input 
impedance,  zin»  when  near  the  earth  as  follows.3 

when  the  antenna  is  horizontal: 

Zin  =  Z,i  +  (-Z21) 

Z21  =  Z|  1  -  Zin  ohmSj  j 


1 


when  the  antenna  is  vertical: 


Zin 

«  Zi  1 

+  (  +  Z2 1 ) 

^2  1 

=  Zin 

-  Zi  1 

ohms^ 

3 

Since  the  antenna  was  pruned  to  first-resonance,  Z\  i  =  Ri 1 ,  and  equations 
2-3  reduce,  respectively,  to: 


Z2  1  = 

Ri  1 

-  Zin 

ohms 

Z2  1  “ 

Zin 

-  R11 

ohms 

The  first-resonant  dipole  self  resistance,  Ri  1  ,  used  in  equations  4  and  5 
can  be  approximated  without  having  to  resort  to  the  computer  NEC  program. 
The  procedure  is  to  use  an  estimated  relative  velocity,  V^Ce),  multiply 
0.5Ao  by  this  V^Ce)  to  get  an  actual  length,  L,  divide  L  by  the  wire 
diameter,  D,  and  solve  equation  4  of  reference  4  for  an  actual  relative 
velocity,  V^(a) .  The  process  is  repeated  (iterated)  until  V^(a)  =  V^(e)  . 

When  each  V^(a)  thus  obtained  is  used  as  the  next  V^(e),  the  solution 
Vf(a)  =  Vf(e)  is  obtained  via  a  relatively  few  iterations  because  of  rapid 
convergence.  Then,  note  L/D  when  V^(a)  =  V^(e),  and  use  this  L/D  in  equa¬ 
tion  9  of  reference  4  to  obtain  Ri  i  .  The  following  is  an  example  solution: 

V^Ce)  =  0.976772  numeric 

0.5Ao  =  2950.718504  inches  at  2.0  MHz 
L  -==  (0.976772)  (2950.718504)  =  2882.179215  inches 
D  =  0.08081  inches  (#12  wire) 

V  (a)  =  1  -  [10.541  logI0(£)  -  4.933]-' 

=  0.976772  =  V  (e)  numeric 


3 


7. ,,  =  R,,  i  73.0  -  (0.252  log,0  (2~ }  +  0.2321"' 

=  72.275  +  j  0.0  ohms 

when  these  L  and  D  values  are  used  in  the  computer  NEC  program,  the 
free-space  solution  at  2  MHz  is 

Z, i  =  72.330  -  j  0.340  ohms 

Thus,  when  L/D  =  35666,  the  above  procedure  gives  a  solution  that  is 
within  0.06  ohms  and  0.27  degrees  of  that  obtained  via  NEC  at  2.0  MHz! 


To  obtain  the  effect  of  frequency  upon  the  results,  the  2-30  MHz  band 
was  divided  into  4  MHz  increments.  Then,  with  the  exception  of  solu¬ 
tions  over  a  perfect  earth,  8  frequencies  are  used  in  this  analysis. 
Where  less  than  8  curves  appear  on  the  imperfect  earth  graphs,  curves 
are  combined  when  the  error  is  less  than  plus-minus  1.5  ohms. 


To  obtain  the  effect  of  earth’s  electrical  properties  upon  the  results, 
5  defined  earths  are  used.  The  arbitrary  properties  are  listed  in 
Table  I. 


TABLE  1 

Earth 

f.  (numeric) 

T  (mhos/meter) 

Poor 

5.0 

0.001 

Typ  ic.al 

13.0 

0.005 

Good 

21.0 

0.010 

Sea  Water 

80.0 

5.000 

Perfect 

1.0 

oo 

The  normalized  height  arguments,  H^,  shown  on  the  graphs  are  f eed  heights 
over  the  earth.  Thus,  the  remote  end  of  one  arm  of  a  vertical  center-fed 
first-resonant  dipole  antenna  is  very  close  to  the  earth  when  =  0.25. 


4 


njjm 


The  curves  the  enclosed  54  figures  were  generated  from  9942  calculated 
data  points.  Thus,  the  enormous  amount  of  required  computer  plus  reduc¬ 
tion  time  precludes  the  inclusion  of  other  antenna  types  in  this  report. 

II.  HORIZONTALLY-POLARIZED  MUTUAL  RESISTANCE. 

The  mutual  resistance,  R2  1  ,  results  are  plotted  on  Figures  2-6,  7-11, 
12-16,  and  17-21  for  height,  H^,  intervals  of  0.001-0.015,  0.015-0.155, 
0.15-0.85,  and  0.85-1.55  wavelengths,  respectively.  Thus,  at  each  height 
interval,  there  are  5  graphs,  one  for  each  defined  earth,  and  the  frequen¬ 
cy  or  frequency  range  is  plotted  on  each  graph. 

With  the  graphs  so  arranged,  some  degree  of  earth  interpolation  is  en¬ 
hanced.  As  an  example,  let  the  earth's  electrical  properties  be  *  10 
and  X  =  0.002  mhos/meter  (between  poor  and  typical  earth).  Using  Figures 
2  and  3  with  H  =  O.OlAo  and  f  =  2.0  MHz,  the  solution  is 
-35.3  <  R2 1  <  +  3.0  ohms.  The  NEC  solution  is  -18.8  ohms. 

These  figures  show  that  ground  (not  sea  water  or  perfect  earth)  mutual 
resistance  is  highly  negative  when  this  antenna  is  near  an  inperfect 
ground  which,  from  equation  4,  increases  the  antenna  resistance  drastic¬ 
ally.  This  conclusion  is  supported  by  field  measurements . 5 ’6 ’  7’  8 

The  results  shown  on  Figure  5  appear  to  be  highly  frequency  sensitive 
when  this  antenna  is  near  sea  water.  The  loss  tangent  of  sea  water 
does  not  exhibit  a  relatively  strong  displacement  current  until  the 
frequency  is  above  1.0  MHz.  The  great  difference  in  R2  1  solutions 
between  26  and  30  MHz  was  not  expected,  and  the  reason  for  this  is 
given  in  the  Summary. 

The  results  shown  on  Figure  6  are  what  one  would  expect.  The  mutual 
resistance,  R2 1  ,  approaches  the  dipole  self  resistance,  Ri  1  ,  when  the 
height,  H^,  approaches  the  dipole  radius,  and  the  dipole  self  resist¬ 
ance,  Ri  1  ,  is  a  function  of  L/D.  For  practical  reasons,  number  12  wire 
(diameter  =  0.08081  inches)  was  used  at  all  frequencies  except  30  MHz 
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where  a  tubing  diameter  of  1.0  inch  was  used.  Equation  4  of  reference  4 
was  used,  as  discussed  in  the  introduction,  to  obtain  first  resonant  I./ll 
ratios  which  are  plotted  on  Figure  52  in  Summa  '  section,  and  shown  on 
Figure  6  at  2,  10,  and  30  MHz. 

When  the  1./0  ratios  shown  on  Figure  6  are  used  in  equation  9  of  reference 
4,  solutions  for  K ■ j  are  0.00  ohms  less  than  NEC  solutions  at  2.0  MHz, 
0.03  ohms  greater  than  NEC  solutions  at  10.0  MHz,  and  0.03  ohms  greater 
than  NEC  solutions  at  30.0  MHz.  The  difference  between  the  NEC  solution 
curves  shown  on  Figure  6  is  a  function  of  dipole  first-resonant  lengths 
of  0.488386\o  at  2.0  MHz.,  0.485980\o  at  10.0  MHz,  and  0.473683Ao  at 
30.0  MHz. 

The  results  plotted  on  Figure  10  show  that  NEC  solutions  are  not  correct 
when  the  dipole  height  is  0.015  <_  ll(  ^  0.03  wavelength  over  sea  water 

and  the  frequency  is  below  2b  MHz.  Surprisingly,  this  height  range  is 
well  within  the  ROM  limit  discussed  in  the  Introduction  leading  to  equa¬ 
tion  1  and  Figure  1.  Since  the  subroutine  gives  reasonable  solutions 
when  other  combinations  of  earth  electrical  properties,  frequency,  and 
height  are  used,  the  error  appears  to  be  philosophical. 

III.  HORIZONTALLY- POLAR I ZED  MUTUAL  REACTANCE . 

The  mutual  reactance,  X; i ,  results  are  ploited  on  Figures  22-2b,  27-31, 
32-30,  and  37-41  for  height,  li ^ ,  intervals  of  0.001-0.015,  0.015-0.155, 
0.15-0.85,  and  0.85-1.55  wavelengths ,  respectively.  Thus,  at  each 
height  interval  there  are  5  graphs,  one  for  each  defined  earth,  and  the 
frequency  or  frequency  range  is  plotted  on  each  graph. 

With  the  graphs  so  arranged,  some  degree  of  earth  interpolation  is  en¬ 
hanced.  As  an  example,  let  the  earth's  electrical  properties  he 
r  =  10  and  t  =  0.002  mhos/meter  (between  poor  and  typical  earth). 

Using  Figures  22  and  23  with  11  =  O.OlAo  and  f  =  2.0  MHz,  the  solution 
is  -90.8  <  X2  1  "s  -  64.6  ohms.  The  NEC  solution  is  -80.7  ohms. 
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The  results  plot  fed  on  Figure  SO  show  ili.it  N  .so  I  ut  ions  are  not  correct 
when  the  dipole  height  is  <1.01  '>  i!  d.ti  t  wave  !  t-uv .1  li  over  sea  water. 

This  is  consistent  with  the  results  obtained  in  So-  t  ion  li,  and  indicate 
that  NEC  -  with  its  existing  espial  ions  -  is  not  vi'id  over  sea  water  at 
HE  when  H .  is  in  this  rev  ion 
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These  figures,  too,  appear  to  be  highly  accurate,  and  vertically- 
polarized  mutual  reactance  is  not  very  significant  at  HF  when 
H ^  '>  0.70  wavelength. 

VI.  SUMMARY. 

While  this  report  is  not  complex,  it  required  so  much  computer  and 
data  reduction  time  that  it  discouraged  any  desire  to  include  solu¬ 
tions  for  other  antenna  lengths.  An  analysis  of  computer  solutions 
indicated  the  validity  range  of  simplified  equations  was  so  narrow 
that  such  an  approach  is  impractical.  This  is  very  apparent  at 
heights,  H-^,  greater  than  0.15  wavelengths  where  solutions  became 
oscillatory.  Hence,  all  of  the  reduced  data  is  plotted  on  the  en¬ 
closed  figures  for  general  use. 

The  accuracy  of  these  results  depends  upon  the  accuracy  of  the  equa¬ 
tions  used  in  program  subroutines.  That  is,  with  the  exception  of 
horizontal  polarization  and  0.01  0.03  wavelength  over  sea 

water,  the  curves  on  the  figures  are  continuous,  and  the  results  are 
highly  predictable.  The  highly  accurate  results  on  Figures  42-51 
indicate  that  computational  errors  will  be  more  related  to  the  Hert¬ 
zian  parallel  electric  component  equation  than  to  the  Hertzian 

perpendicular  electric  component  equation  used  in  the  program 

Sommerfeld  subroutine. 

As  noted  above,  solutions  plotted  on  Figures  12-21  and  32-41  are 
oscillatory  when  >  0.15  wavelength.  The  individual  figures  suggest 
that  R2 1  and  X2  1  solutions  are  highly  sensitive  to  frequency  when  the 
earth's  permittivity  is  highly  conductive.  When  the  figures  are  re¬ 
viewed  collectively,  frequency  is  less  important  than  the  earth's  per¬ 
mittivity.  This  leads  to  the  conclusion  that  the  frequency,  per  se, 
effect  is  more  one  of  length,  L, effect  where  length  at  first  resonance 
is  a  function  of  the  length-to-diameter ,  L/D,  ratio  equation  4  in  Ref¬ 
erence  4. 
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The  practical  L/D  ratios  used  in  this  report  are  plotted  on  Figure  52, 
where  No.  12  wire  (D  =  0.08081  inches)  was  used  at  2-26  MHz  and  1.0 
inch  diameter  tubing  was  used  at  30  MHz.  Thus,  the  6  MHz  first-reson¬ 
ant  dipole  is  0.3%  shorter  than  the  2  MHz  first-resonant  dipole,  the 
10  MHz  first-resonant  dipole  is  0.5%  shorter  than  the  2  MHz  first-re- 
sonant  dipole,  the  14  MHz  first-resonant  dipole  is  0.6%  shorter  than 
the  2  MHz  first-resonant  dipole,  the  18  MHz  first-resonant  dipole  is 
0.7%  shorter  than  the  2  MHz  first-resonant  dipole,  the  22  MHz  first- 
resonant  dipole  is  0.8%  shorter  than  the  2  MHz  first-resonant  dipole, 
the  26  MHz  first-resonant  dipole  is  0.9%  shorter  than  the  2  MHz  first- 
resonant  dipole,  and  the  30  MHz  first-resonant  dipole  is  3 . 1%  shorter 
than  the  2  MHz  first-resonant  dipole.  This  behavior  is  apparent  near 
the  max-min  regions  on  Figures  12-21  and  32-41.  It  also  gives  an 
explanation  for  the  relatively  large  capacitive  solutions  at  30  MHz 
on  Figures  2-5  when  the  first-resonant  dipole  is  very  near  the  ground. 

If  the  30  MHz  first-resonant  dipole  had  been  made  of  No.  12  wire,  it 
would  have  been  2.1%  longer  than  that  of  the  2  MHz  first-resonant 
dipole,  the  NEC  solution  on  Figure  4  would  have  been  -82.9  ohms  when 
=  0.002  wavelength,  and  the  NEC  solution  on  Figure  24  would  have 
been  -210.6  ohms  when  H^  =  0.002  wavelength  (approximate  the  26  MHz 
solutions )  On  the  other  hand,  if  the  30  MHz  first  resonant  dipole 
L/D  ratio  had  been  the  same  as  that  of  the  2  MHz  first  resonant  dipole, 
the  NEC  solution  on  Figure  4  would  have  been  -69.9  ohms  when  H^  =  0.002 
wavelength,  and  the  NEC  solution  on  Figure  24  would  have  been  -202.8 
ohms  when  H^  =  0.002  wavelength  (approximate  the  14  MHz  solutions I) 
Therefore,  at  low  antenna  heights,  solutions  are  highly  dependent 
upon  dipole  length. 

The  behavior  of  mutual  R2  1  and  X2 1  vs.  H^  for  the  5  defined  near  earths 
can  be  presented  in  general  terms  when  frequency  {  and  L/D]  dependence 
is  eliminated.  The  horizontally-polarized  results  at  all  8  frequencies 
were  averaged  for  each  earth,  and  the  average  results  vs.  H^  are  plotted 
on  Figures  53  and  54.  These  results  show: 
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1,  The  inaccuracy  of  NEC  in  the  0.01  £  _<  0.03  region  over 

sea  water . 

2.  Both  R21  and  X21  have  maximum  values  as  a  function  of  earth's 
electrical  properties , 

3,  The  height,  H^,  at  which  maximum  Rn  and  X2 1  ,  occur  depends 
inversely  upon  how  good  the  earth  is  as  a  conductor. 

4.  With  the  exception  of  sea  water  and  perfect  earth,  the  mutual 
impedance  terras,  R2 1  and  Xj  1  ,  are  highly  negative  when  the  horizontal 
first-resonant  dipole  is  very  near  the  earth. 

While  the  results  presented  in  this  report  are  not  precise,  they  are  as 
accurate  as  one  can  expect  without  resorting  to  solutions  involving  a 
large  number  of  L/D  ratios.  The  problem  is,  as  it  turns  out,  that 
solutions  ate  as  much  a  function  of  first-resonant  dipole  length  as 
they  are  to  the  electrical  properties  of  the  earth  beneath  the  dipoles. 

I  am  indebted  to  a  number  of  PED  personnel.  Mr.  Danny  Fink  set  up  the 
programs.  Ms.  Lee  Ann  Sampson  and  SP5  Virgil  Brown  were  the  terminal 
operators.  SSG  Robert  Pulliam,  SP5  Alvin  Mack,  and  Mr.  Steve  Aubrey 
collated  the  stacks  of  printouts.  Without  the  coordinated  effort,  of 
all,  this  lrngthy  report  would  have  been  impossible. 
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